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Abstract. Using the United Nations Commodity Trade Statistics Database we construct the Google 
matrix of the world trade network and analyze its properties for various trade commodities for all countries 
and all available years from 1962 to 2009. The trade flows on this network are classified with the help of 
PageRank and CheiRank algorithms developed for the World Wide Web and other large scale directed 
networks. For the world trade this ranking treats all countries on equal democratic grounds independent 
of country richness. Still this method puts at the top a group of industrially developed countries for trade 
in all commodities. Our study establishes the existence of two solid state like domains of rich and poor 
countries which remain stable in time, while the majority of countries are shown to be in a gas like phase 
with strong rank fluctuations. A simple random matrix model provides a good description of statistical 
distribution of countries in two-dimensional rank plane. The comparison with usual ranking by export and 
import highlights new features and possibilities of our approach. 

PACS. XX.XX.XX No PACS code given 



1 Introduction 

The analysis and understanding of world trade is of pri- 
mary importance for modern international economics |2 . 
Usually the world trade ranking of countries is done ac- 
cording their export and/or import counted in USD [3]. 
In such an approach the rich countries naturally go at the 
top of the listing simply due the fact that they are rich 
and not necessary due to the fact that their trade network 
is efficient, broad and competitive. In fact the trade be- 
tween countries represents a directed network and hence 
it is natural to apply modern methods of directed net- 
works to analyze the properties of this network. Indeed, 
on a scale of last decade the modern society developed 
enormously large directed networks which started to play 
a very important role. Among them we can list the World 
Wide Web (WWW), Facebook, Wikipedia and many oth- 
ers. The information retrieval and ranking of such large 
networks became a formidable challenge of modern soci- 
ety. 

An efficient approach to solution of this problem was 
proposed in [3| on the basis of construction of the Google 
matrix of the network and ranking all its nodes with the 
help of the PageRank algorithm (see detailed description 
in [5]). The elements Gij of the Google matrix of a network 
with N nodes are defined as 



= aS,j + (1 - a)/N 



(1) 



where the matrix S is obtained by normalizing to unity 
all columns of the adjacency matrix Aij, and replacing 
columns with only zero elements by 1/iV. Usually for the 



WWW an element Aij of the adjacency matrix is equal 
to unity if a node j points to node i and zero otherwise. 
The damping parameter a in the WWW context describes 
the probability (1 — a) to jump to any node for a random 
surfer. The value a — 0.85 gives a good classification for 
WWW [5]. By construction the Google matrix belongs 
to the class of Perron-Frobenius operators and Markov 
chains [5], its largest eigenvalue is A = 1 and other eigen- 
values have |A| < a. According to the Perron-Frobenius 
theorem the right eigenvector, called the PageRank vec- 
tor, has maximal A = 1 and non-negative elements that 
have a meaning of probability P{i) attributed to node 
i. Thus all nodes can be ordered in a decreasing order of 
probability P{i) with the corresponding increasing PageR- 
ank index K{i). The presence of gap between A = 1 and 
|A| = a ensures a convergence of a random initial vec- 
tor to the PageRank after about 50 multiplications by 
matrix G. Such a ranking based on the PageRank algo- 
rithm forms the basis of the Google search engine [S] . It 
is established that a dependence of PageRank probabil- 
ity P{i) on rank K{i) is well described by a power law 
P{K) oc l/iC^'" with /3„ « 0.9. This is consistent with 
the relation /3i„ = ^/{^lin — 1) corresponding to the aver- 
age proportionality of PageRank probability P{i) to its in- 
degree distribution Win{k) oc l/A:^*" where k{i) is a num- 
ber of ingoing links for a node i [6, 5J . For the WWW it is 
found that for the ingoing links w 2.1 (with /3i„ « 0.9) 
while for out-degree distribution Wout of outgoing links a 
power law has the exponent /Xoiit ~ 2.7 [ZUS]. We note that 
PageRank is used for ranking in various directed networks 
including citation network of Physical Review [^[TU] and 
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for rating of the total importance of scientific journals [11] . 




5q'" q'" 

Fig. 1. (Color online) Image of money mass matrix M (top), 
Google matrix G (middle) and inverse Google matrix G* (bot- 
tom) for all commodities (left column) and crude petroleum 
(right column) for year 2008 with all world countries = 227 
from the UN COMTRADE [T]. Matrix elements g, for M^j, 
Gi,j or G*j, are shown by color changing from to a corre- 
sponding maximum value g-max- AH three matrices are shown 
in the basis of PageRank index K (and K') of matrix G, respec- 
tively for all commodities (left) and crude petroleum (right), 
which correspond to x, y-axis with 1 < K, K' < N . Here we 
use a — 0.5 for matrix G and its PageRank index K and the 
same a for G*. 



The PageRank performs ranking determined by in- 
going Hnks putting at the top most known and popular 
nodes. However, in certain networks outgoing links also 
play an important role. Recently, on an example of proce- 
dure call network of Linux Kernel software, it was shown 
[T^ that a relevant additional ranking is obtained by tak- 
ing the network with inverse link directions in the ad- 



jacency matrix corresponding to Aij — !■ — Aji and 
constructing from it an additional Google matrix G* ac- 
cording to relation ^ at the same a. The examples of 
matrices G and G* for the world trade network are shown 
in Fig. [TJ The eigenvector of G* with eigenvalue A = 1 
gives then a new inverse PageRank P*{i) with ranking 
index K*{i). This ranking was named CheiRank 13] to 
mark that it allows to chercher I 'information in a new way. 
While the PageRank rates the network nodes in average 
proportionally to a number of ingoing links, the CheiRank 
rates nodes in average proportionally to a number of out- 
going links. The results obtained in [T^[T5] confirm this 
proportionality with the exponent /3out = l/(7o«t ~ !)• 

Since each node belongs both to CheiRank and PageR- 
ank vectors the ranking of information flow on a directed 
network becomes two-dimensional. While PageRank high- 
lights how popular and known is a given node, CheiRank 
highlights its communication and connectivity abilities. 
The examples of Linux and Wikipedia networks show that 
the rating of nodes based on PageRank and CheiRank al- 
lows to perform information retrieval and to characterize 
network properties in a qualitatively new way jl2lll3| . 

In this work we apply CheiRank and PageRank ap- 
proach to the World Trade Network (WTN) using the 
enormous and detailed United Nations Commodity Trade 
Statistics Database (UN COMTRADE) [J. Using these 
data we analyze the world trade flows both in import and 
export for all commodities for all years 1962 - 2009 avail- 
able there at SITCl and HS96 databases. We also per- 
formed analysis for specific commodities taken from SITC 
Rev. 1 database, mainly for year 2008: crude petroleum 
(Sl-33101, "Crude petroleum"), natural gas (Sl-3411, "Gas, 
natural"), barley (Sl-0430, "Barley, unmilled"), cars (Sl- 
7321, "Passenger motor cars, other than buses"), food (Sl- 
0, "Food and live animals"), cereals (Sl-04, "Cereals and 
cereal preparations" ) . Their codes and official UN names 
are given in brackets. In few cases, when certain coun- 
tries were non-reporting their export, we complemented 
the WTN data from the import database. 

For a given year we extract from the UN COMTRADE 
money transfer (in USD) from country j to country i that 
gives us money matrix elements Mij (for all types of com- 
modities noted above). These elements can be viewed as 
a money mass transfer from j to i. In contrast to the ad- 
jacency matrix Aij of WWW, where all elements are only 
or 1, here we have the case of weighted elements. This 
corresponds to a case when there are in principle multiple 
number of links from j to i and this number is propor- 
tional to USD amount transfer. Such a situation appears 
for rating of scientific journals [TT], Linux PCN [12^ and for 
Wikipedia English articles hyperlink network 13 , where 
generally there are few citations (links) from a given ar- 
ticle to another one. In this case still the Google matrix 
is constructed according to the usual rules and relation 
^ with Sij = Mij/nij and Sij ~ ^/N, if for a given 
j all elements Mij = 0. Here nij = J^i ^ij is the total 
export mass for country j. The matrix G* is constructed 
from transposed money matrix with Sij ~ ^Iji/ Tli-^ji- 
In this way we obtain the Google matrices G and G* of 
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WTN which allow to treat all countries on equal grounds 
independently of the fact if a given country is rich or poor. 
A similar choice was used in rating of scientific journals 
[IT] . PCN Linux [12] and Wikipedia network [13]. The 
main difference appearing for WTN is a very large varia- 
tion of mass matrix elements My related to the fact that 
there is very strong variation of richness of various coun- 
tries. Due to these reason we think that it is important to 
use the ranking based on the Google matrix which treats 
in a democratic way all world countries that corresponds 
to the democratic standards of the UN. For the WTN 
CheiRank and PageRank are naturally linked to export 
and import flows for a given country and hence it is very 
natural to use these ranks for characterization of country 
trade abilities. The Google matrix can be constructed in 
the same way not only for all commodities but also for a 
given specific commodity. 

We note that recently the interest to the analysis of the 
world trade as a network becomes more and more pro- 
nounced with a few publications in this area [H1IT51ITB1 
[TT] . Thus, the global network characteristics were consid- 
ered in P^[T5] , degree centrality measures were analyzed 
in [IB] and time evolution of network global characteristics 
was studied in [17]. Topological and clustering properties 
of multinetwork of various commodities were discussed 
in [IS]. Here we present a systematic study of directed 
WTN on the basis of new combination of PageRank and 
CheiRank methods using the Google matrix constructed 
for the enormous UN COMTRADE database. 

The paper is composed as follows: in Section 2 we 
describe the global properties of the Google matrix of 
WTN, in Section 3 we analyze distribution of countries in 
PageRank-CheiRank plane for all time period 1962 - 2009 
and propose a random matrix model of WTN (RMWTN) 
which describes the statistical properties of this distribu- 
tion in the case of all commodities; comparison with rank- 
ing based on import and export for various commodities is 
presented in Section 4; discussion of the results is given in 
Section 5. More detailed information and data are given 
in Appendix and at the website T^. 

2 Properties of Google matrix of WTN 

An example of the Google matrix of WTN in 2008 is shown 
in Fig. [T] for all commodities and crude petroleum. The 
matrices G and G* are shown in the bases where all coun- 
tries are ordered by the PageRank index K of matrix G 
constructed for corresponding commodity (left and right 
columns). The matrix elements of G are distributed over 
all N values being roughly homogeneously in K, even if 
the left top corner at small K, K' values is filled in a more 
dense way. In contrast the density drops at large values 
of K' . Such a structure is visible both for all commodi- 
ties and crude petroleum but clearly the global density is 
smaller in the later case since there are less number of 
links there (see data in next Section). The structure of 
G* is approximately the same (we will see in next Section 
that rich countries are located at low K, K' values). In 
contrast to G and G* the structure of money matrix M is 



rather different. For all commodities matrix elements drop 
very rapidly at large values of K and K' that corresponds 
to the fact that the main amount of world money circu- 
lates only between rich countries with top ranks K. In 
contrast to that for crude petroleum the matrix elements 
of M are located at intermediate K values. Indeed, in this 
case PageRank index K orders countries by their crude 
petroleum trade where richest countries are not necessar- 
ily at the top ranks. 




Fig. 2. (Color online) Probability distributions of PageR- 
ank P(ii:), CheiRank P'{K'), ImportRank P(k), and Ex- 
portRank P*{K*) are shown as function of their indexes in 
logarithmic scale for all commodities (top panel) and crude 
petroleum (bottom panel) for WTN in 2008 with N = 227. 
Here P(K) and P*{K*) are shown by red and blue curves re- 
spectively, for a = 0.5 (solid curves) and a = 0.85 (dotted 
curves); P{K) and P*{K*) are displayed by dashed red and 
blue curves respectively. For both commodities the distribu- 
tions P{K) and P*{K*) follow a power law dependence like 
P oc l/K^ (see text), the Zipf law is shown by the straight 
dashed line with /J = 1 in top panel. 



From the Google matrices G and G* we find the proba- 
bility distributions PageRank P{K) and CheiRank P*{K*) 
which are shown in Fig. [2] for the same commodities as 
Fig. [TJ One of the main features of these distributions is 
that both P{K) and P*{K*) depend on their indexes in a 
rather similar way from, that is in contrast to the results 
found for the WWW [HIS], PCN Linux [HI and Wikipedia 
network [131, where these distributions are different hav- 
ing different exponents /? in the power law decay. Here, 
up to fluctuations, we have /3i„ = 13 out = P- The size of 
WTN is rather small compared to usual sizes of WWW, 
Linux or Wikipedia networks. However, still we find that 
the power law gives a quite good fit of our data. The fit 
gives (3 = 1.17 ± 0.015 at a = 0.85 and (3 = 0.63 ± 0.01 
at a = 0.5 (for all commodities) and (3 = 0.92 ± 0.02 
and (3 = 0.51 ±0.01 respectively (for crude petroleum) for 
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all 227 countries in Fig. [21 For the fit of top 100 coun- 
tries we have respectively /3 = 1.15 ± 0.03 (a = 0.85) 
and /3 = 0.75 ± 0.008 (a — 0.5) for all commodities and 
(3 = 1.22±0.015 {a = 0.85) and /3 = 0.70±0.008 {a = 0.5) 
for crude petroleum. There is a certain change of the ex- 
ponent with a decrease of the fit interval which, however, 
is not very large. We attribute this to visible deviations at 
the tail of K, K* with small countries (see discussion in 
next Section). In average the exponent value is not very 
far from the value /3 = 1 corresponding to the Zipf law 
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Fig. 3. (Color online) Spectrum of the eigenvalues A of the 
Google matrix at a = 1 in complex plane for WTN in 2008 
with N = 227 countries for all commodities, food, cereals and 
barley (from top to bottom); all eigenvalues are shown for each 
commodity; unit radius circle is also shown. Only the bottom 
case have quasi-degenerate eigenvalues close to the circle with 
3 values A = 1., 0.99987, 0.991 and two values close to A = — 1; 
other cases have a significant gap separation from |A| = 1. 



It is useful to compare the behavior of probabilities 
P{K) and P*{K*) with respective ranking related to im- 
port and export. To do that we rank the countries by prob- 
ability import P{K) defined as a ratio of import in USD 
for a given country K to the total world import in USD 
for a given year with ordering of all countries in decreas- 
ing probability order index of ImportRank K. By con- 
struction we have J2k^i^) — ^ ^^'^ P{^) = "^if A^t, 
where m^^ is the import mass of a given country K and 
Mt = j ^ij is the total world money mass for a given 
year. In the same way we construct export probability 
P*{K*) with the ExportRank K*. The dependence of 
these probabilities on their indexes is shown in Fig. [21 In 
the range of 1 < K, K* < 50 it can be well described by a 
power law with /3 = 1.01 ± 0.03 for all commodities corre- 
sponding to the Zips law (for crude petroleum we obtain 
for this range (3 = 1.43 ± 0.07). At larger values of or- 
der index we find a sharp drop with an exponential type 
decay on the tail. For crude petroleum this exponential 
decay starts at smaller values of K due to a significantly 



smaller total number of links that gives an increase of /3 for 
the range of 50 countries. The exponential decay at large 
K results from a strong variation of richness of countries 
which changes more than by four orders of magnitude. 
From the comparison of ranks shown in Fig. [51 it is clear 
that PageRank and CheiRank give more equilibrated and 
democratic description of trade flows. 

We should note that due to a small size of the WTN 
the fluctuations are stronger compared to large size net- 
works like the WWW. It is especially visible for specific 
commodities where the total number of links is by fac- 
tor 30 smaller then for all commodities (see next Sec- 
tion). These fluctuations are smaller for the damping fac- 
tor value a = 0.5 in agreement with the results presented 
in pni^ . In fact this a value was also used in [5] for 
PhysRev citation network. Due to that reasons in next 
Sections we show data for ranking at a — 0.5. 

Finally let us discuss the spectrum A of the Google ma- 
trix which follows from the equation for right eigenvectors 

1pm{i)- 

^GijIprnU) ^ ^milmii) ■ (2) 
j 

It is known that the dependence on a is rather simple: all 
eigenvalues, except one with A = 1, are multiplied by a 
[5]. Due to that we show the spectrum of G at a = 1 in 
Fig. [21 Compared to the spectrum studied for other net- 
works (see examples in |121I221I231I24II251I26| 1 we find that 
the WTN spectrum is very close to real line especially 
for three top commodities in Fig. [31 We explain this by 
the fact that here an average number of links per coun- 
try is very large for these commodities and that the ma- 
trix elements are not very far from the symmetric relation 
Mij = Mji at which the spectrum is real. We only note 
that for barley the spectrum has quasi-degeneracy at A = 1 
that signifies the existence of slow relaxation modes. We 
attribute this to the fact that there are certain countries 
which practically do not use barley that leads to appearing 
of isolated subspaces with corresponding quasi-degenerate 
modes. We will return to the discussion of spectrum prop- 
erties of G in next Section. 



3 CheiRank versus PageRank for WTN 

We start from examples of distributions of countries in the 
PageRank-CheiRank plane shown in Fig. [4l for 5 different 
commodities in 2008. The first case of all commodities cor- 
responds to trade flows between countries integrated over 
all type of products. Even if the Google matrix approach 
is based on a democratic ranking of international trade, 
being independent of total amount of export-import for a 
given country, we still find at the top ranks K and K* the 
group of industrially developed countries (see more details 
in Table 1 in Appendix). This means that these countries 
have efhcient trade networks with optimally distributed 
trade flows. Another pronounced feature of global distri- 
bution is that it is concentrated along the main diagonal 
K = K* . This feature is not present in other networks 
studied before (e.g. PCN Linux [T21 and Wikipedia [13]). 
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The origin of this density concentration is based on sim- 
ple economy reason: for each country the total import is 
approximately equal to export since each country should 
keep in average an economic balance. Thus for a given 
country its trade is doing well if its K* < K so that the 
country exports more than it imports. The opposite rela- 
tion K* > K corresponds to a bad trade situation. We also 
can say that local minima in the curve of (K* — K) vs. K 
correspond to a successful trade while maxima mark bad 
traders. In 2008 most successful were China, Rep. of Ko- 
rea, Russia, Singapore, Brasil, South Africa, Venezuela (in 
order of K ior K < 50) while among bad traders we note 
UK, Spain, Nigeria, Poland, Czech Rep., Greece, Sudan 
with especially strong export drop for two last cases. The 
comparison of our ranking with the import-export ranking 
will be analyzed in next Section. 

Even if there is a concentration of density along the 
main diagonal (Fig.|4^) we still have a significant broaden- 
ing of distribution especially at middle values of if '-^ 100. 
This means that the gravity model of trade, often used in 
economy (see e.g. pill6|). has only approximate validity. 
Indeed, in this model the mass matrix Mij is symmetric 
that would placed all countries on diagonal K = K* that 
is definitely not the case. 

If we now turn to the distribution of countries for 
a trade in a specific commodity than it becomes abso- 
lutely clear that the symmetry approximately visible for 
all commodities is absolutely absent: the points are scat- 
tered practically over the whole square N xN. The reason 
of such a strong scattering is clear: e.g. for crude petroleum 
some countries export this product while other countries 
import it. Even if there is some flow from exporters to 
exporters it remains relatively low (see more discussion in 
next Section). This makes the Google matrix to be very 
asymmetric. Indeed, the asymmetry of trade flow is well 
visible in Fig. Eh where arrows show the trade directions 
between countries within top 40 x 40 ranks for barley. 

It is also useful to use 2DRank K2 discussed in [T3], 
which orders all nodes according to the order of their ap- 
pearance inside squares of size K x K going from K ~ 
1,2,3, ... to N for each of four specific commodities shown 
in Fig. m In a certain sense top countries in 2DRank K2 
are those which are active traders even not being among 
large exporters or importers of this product (all ranks 
for commodities of Fig. 2] are given in Tables 1-5 in Ap- 
pendix). As an example, we note Singapore which is at 
the third position in K2 (Table 2): it is a small country 
which cannot export or import a large amount of the com- 
modity, but its trade network is very active redistributing 
flows between various countries that places it at a high K2 
rank. 

The images of Fig. fallow to understand qualitatively 
the reasons of density concentration around diagonal K — 
K* for the case of all commodities: this trade is composed 
from hundreds of specific commodities which behave ran- 
domly and the averaging over them gives effective coarse- 
graining and produces a certain symmetry for matrix ele- 
ments due to the central limit theorem for a sum of many 
positive contributions. The fact that the increase of coarse- 
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Fig. 4. (Color online) Country positions in PageRank- 
CheiRank plane {K,K*) for world trade in various commodi- 
ties in 2008. Each country is shown by circle with its own flag 
(for a better visibility the circle center is slightly displaced 
from its integer position {K,K*) along direction angle ■7r/4). 
The panels show the ranking for trade in the following com- 
modities: all commodities {a,b); crude petroleum {c,d); natural 
gas (e, /); barley (g, h); cars {i,j). Left column shows a global 
scale with all 227 countries, while right column gives a zoom 
in the region of 40 x 40 top ranks. For barley in panel (h) the 
links between countries inside the selected region are shown by 
arrows. 



graining cell gives more and more symmetry is well seen in 
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Fig. 5. (Color online) Time evolution of the number of coun- 
tries (TV, full black curve), average number of links per country 
{{Nl}) for all commodities (dashed curve) and crude petroleum 
(points for five years), total amount of money (Mt = ^ Mij, 
red curve). The scale of A'' and (Nl) is shown on left side, while 
Mt values, in $ USD, are given in logarithmic scale on the right 
side. 



Fig. 12] where the spectrum becomes more and more close 
to a real one, and hence there is more and more symmetry 
in elements Gij, when we go from barley to cereals, food 
and all commodities. 

We will return to the analysis of specific country rank- 
ing in the next Section while now we turn to analysis of 
time evolution of WTN. 

The variation of global parameters of WTN during the 
database period 1962 - 2009 is shown in Fig. [5] The num- 
ber of countries is increased by 38%, while the number of 
links per country for all commodities is increased in to- 
tal by 140% with a significant increase from 50% to 140% 
during the period 1993 - 2009 corresponding to economy 
globalization. At the same time for a specific commod- 
ity the average number of links per country remains on a 
level of 3-5 links being by a factor 30 smaller compared to 
all commodities trade. During the whole period the total 
amount Mt of trade in USD shows an average exponential 
growth by 2 orders of magnitude. 

To understand the physical properties of the WTN we 
consider the distribution of money mass transfer matrix 
elements Mij shown versus their transposed values Mji in 
Fig- EI This distribution is symmetric by the construction. 
In the case of symmetric matrix Mij , corresponding to the 
gravity model of trade or undirected network, all elements 
should be located on one diagonal line that is definitely 
not the case. For crude petroleum the distribution is even 
more broad showing definite absence of symmetry of Mij . 
In fact for all commodities the distribution forms a rather 
broad cone which form remains stable in time according to 
the comparison of data in 1962 and 2008 years (the den- 
sity is higher in the later case since there are more coun- 
tries). Keeping in mind that according to data of Fig. [5] 
we have the Zipf law for P{K) we propose the random 
matrix model of WTN (RMWTN) with the mass matrix 
elements given by 

Mij = e.tej/ij , (3) 




M. . M. . 



Fig. 6. Money mass transfer matrix elements Mij are shown 
versus their transposed values Mj,i for all commodities of WTN 
in 1962 (top left panel) and 2008 (top right panel). Bottom 
left panel shows the matrix elements for crude petroleum of 
WTN in 2008; bottom right panel shows the same quantities 
for random matrix model of WTN. Four panels show 5 orders of 
magnitude in logarithmic scales starting from maximum values 
of Mij . In the case of WTN (top and bottom left panels) matrix 
elements are taken from the UN COMTRADE database and 
are expressed in USD, right bottom panel is built from one 
random realization with Mij = eiej/ij (see text). Here A'^ = 
164 for 1962 data; iV = 227 for 2008 data and RMWTN. 



where are random numbers homogeneously distributed 
in [0, 1] interval and i,j are indexes in the ImportRank in- 
dex K. The distribution given by this simple model repro- 
duces quite well the actual distribution found for all com- 
modities (see right panels in Fig. [5]). With this RMWTN 
distribution of Mij we construct the Google matrices G 
and G* according to the usual recipes ([TJ and then deter- 
mine the distribution of points in {K,K*) plane. 

To have a statistical comparison between the RMWTN 
and real WTN data we construct the density distribution 
of countries in the plane {K* — K, K* + K) using all avail- 
able years 1962 - 2009 at the UN COMTRADE database 
for all commodities. The coarse-grained distribution of 
about 10^ WTN data points is shown in Fig.jT] We present 
the data directly in {K*—K, K*-\-K) plane (top left panel) 
and in rescaled variables {{K*-K)/N, {K*+K)/N) plane, 
which takes into account that the number of countries 
grown by 38% during this time period. The distribution 
has a form of spindle with maximum density at the ver- 
tical axis K* — K = Q. We remind that good exporters 
are on the left side of this axis at K* — K < Q, while the 
good importers (bad exporters) are on the right side at 
K* -K>Q. 

The comparison of WTN data with the results pro- 
duced by RMWTN model ^ are shown in bottom panels 
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Fig. 7. (Color online) Spindle distribution for WTN of all 
commodities for all countries in the period 1962 - 2009 shown 
in the plane of {K* — K, K* + K) (left top panel, coarse- 
grained data in 3 X 3 cell size) and in the rescaled plane 
{(K* — K)/N, {K* + K)/N) (right top panel, coarse-graining 
inside each of 76 x 152 cells, which is approximately the same 
number as in top left panel); data from the UN COMTRADE 
database. Bottom left panel: zoom of top right panel; bottom 
right panel: data from 100 realisations of RMWTN model © 
with iV = 227 as for WTN size in 2008. 



of Fig. [71 there is a good agreement between both without 
any fit parameters for the half of all countries with top 
ranks [K + K* < N). For countries with K + K* > N 
RMWTN model does not succeed to describe correctly the 
upper part of spindle distribution found for the WTN and 
hence further improvements of the RMWTN are needed. 
However, a simple description of the distribution for a half 
top countries is rather successful. 

A remarkable feature of the WTN spindle distribu- 
tion of Fig. [7] (top right) is the appearance of high den- 
sity domains sA K* — K Q with K -\- K* sa 1 and 
K+K* K, 2N. They give an impression of two solid phases 
emerging in these two regions while the other part looks 
like a gas phase. This view gets additional confirmation 
by data of Fig. |S] where we present the velocity square 
(Z\w)^, averaged over the whole period 1962 - 2009, as 
a function oi K + K* . This local quantity is defined as 
(Z\«)2 = {K{tl) - K{t - + {K*(t) - K*{t - via 
a one year displacement of a given country in (K^K*) 
plane with further averaging over all times and all coun- 
tries. These data clearly show that for K + K* < 20 we 
have very small square velocity (small effective temper- 
ature) corresponding to a solid phase of rich countries. 




Fig. 8. (Color online) Top panel shows velocity square Av^ 
as a function of K + K* for all countries and all years [all 
commodities data). Gray circles represent all values of Av^ , 
red curve shows the value of Av^ averaged over cases with 
fixed K -\- K* , blue curve shows the average of the red curve 
data in the interval [K + K* - 10, K + K* + 10]. In the bottom 
panel the number of cases p{K + K*) at a given K + K* is 
shown as a function of K + K* . 



while for K + K* > 20 we have large square velocity (large 
effective temperature) corresponding to a gas phase with 
rapid rank fluctuations. There is a similar visible drop of 
temperature at another limit of most poor countries with 
K + K* w 2A^ which indicates a formation of solid phase 
of poor countries (the data are not so exact for this region 
due to variation of number of UN countries with time). 

The presence of solid phase of rich countries and gas 
phase of other countries is also visible from analysis of 
rank variation in time for individual countries shown in 
Fig. O for K, K* < 10 the curves are almost flat while for 
K, K* > 10 we see strong fluctuation of curves. It is in- 
teresting to note that sharp increases in K mark crises in 
1991, 1998 for Russia and in 2001 for Argentina (import 
is reduced in period of crises). We also see that in recent 
years the solid phase is perturbed by entrance of new coun- 
tries like China and India. However, the results presented 
in Fig. [To] for the variation of square velocity with time for 
three regions oi K + K* show that the top 10, and even 
top 20, countries have rather small velocities Av"^, com- 
pared to those with {K + K*)/2 K > 20. For K < 20 
we have Av^ which remains constant in time. In a certain 
sense it looks that the countries with 20 < iiT < 40 pro- 
tect those with 1 < K < 20 (approximately corresponding 
to G-20 major economies [27 ), so that their temperature 
at 1 < iiT < 20 remains unaffected even by a very larger 
fluctuation well visible for the range 81 < K + K* < 120 
during the period of 1992 - 1998 with a few financial crises 
of Black Wednesday, Mexico crisis, Asian crisis and Rus- 
sian crisis. 

The presented results for distribution of countries and 
analysis of their time evolution in the PageRank-CheiRank 
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year year 

Fig. 9. (Color online) Time evolution of CheiRank and PageR- 
ank indexes K, K* for some selected countries for all commodi- 
ties. The countries shown in top panels are: Japan (jp- black), 
France (fr-red), Fed. Rep. of Germany and Germany (de - both 
in blue), Great Britain (gb - green), USA (us - orange) [curves 
from top to bottom in 1962 in left top panel]. The countries 
shown in bottom panels are: Argentina (ar - violet), India (in 
- dark green), China (cn - cyan), USSR and Russian Fed. (ru - 
both in gray) [curves from top to bottom in 1975 in left bottom 
panel] . 



600 1 




°lWo 1970 1980 1990 2000 2010 

years 

Fig. 10. (Color online) Time evolution of velocity square Av^ 
for all commodities averaged over five years interval. In ad- 
dition Av'^ is averaged over countries in the following inter- 
vals: 1 < K + K* < m (blue curve), 41 < K + K* < 
(red curve), SI < K -\- K* < 120 (black curve) in top panel; 
1<K + K'' <20 (blue curve), 21 < K + K* < 40 (red curve), 
41 < K + K* < 60 (black curve) in bottom panel. 



plane confirm a weW knovifn statement that "the poor stay 
poor and the rich stay rich". 

Finally let us discuss an additional parameter which 
characterizes the correlation between PageRank and Chei- 
Rank vectors. The correlator between PageRank and Chei- 
Rank is defined as 

K^NY,PiKi^))P*{K*i^))-^, (4) 

i 

and in a similar way the correlator between ImportRank 
and ExportRank is given by 

k^Nj2P{Kiz))P*iK*{i))-l. (5) 

i 

Recently it has been found that there are networks with 
small correlator, like PCN Linux [T^, and large correlator, 
as Wikipedia [13]. Here we find that for all commodities 
we have large values of k and k, which have rather sim- 
ilar dependence on time (see Fig. [TT|) . In contrast, there 
are almost zero or even negative correlations for crude 
petroleum. Indeed, for crude petroleum there is no correla- 
tion between export and import while for all commodities 
they are strongly correlated. 



4 Comparison with Import - Export ranking 

It is important to compare rating based on PageRank and 
CheiRank with the useful way of country rating based on 
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years 

Fig. 11. (Color online) Time evolution of correlators of 
PageRank-CheiRank (k) and ImportRank-ExportRank (k). 
All commodities are shown by solid red curve for k, and solid 
black curve and dashed blue curve for k with a — 0.5 and 
a = 0.85 respectively. Correlators for crude petroleum with 10 
years of separation are shown in red squares for R and black 
circles for k. 



ImportRank and ExportRank (see e.g [3]). With this aim 
we present the distribution of country positions on the 
planes {K,K) and {K*,K*) shown for top 100 for the 
same commodities as in Fig. |3]for year 2008. For all com- 
modities there is a clear correlation between PageRank 
and ImportRank since the distribution of points is cen- 
tered along the diagonal K = K. It starts to spread only 



L.Ermann, D.L.Shepelyansky: Google matrix of the world trade network 



9 



around K K ^ 30. At the same time for CheiRank 
and ExportRank such a spreading from diagonal starts 
significantly earlier at K* w K* « 10. 

For other commodities shown in Fig. [T2] the correla- 
tions between ranking based on Google matrix and cor- 
responding Export or/and Import ranking are practically 
absent showing very broad scattering of points around al- 
most the whole plane. Only for cars there is a certain level 
of correlation for approximately the first 10 ranks. Natural 
products like crude petroleum, natural gas and agriculture 
products like barley show no correlations. 

The similar conclusions can be also drawn from the 
comparison of country distributions in the plane {K,K*) 
(Fig. El) and in the plane {K,K*) (Fig. [TS]), which show 
data on the same scales. Clearly, the distributions are 
rather different and only for all commodities we can see 
visible correlations (we note that appearance of ordered 
short line segments in panels (c,g) around K Ki K 200 
is due to degeneracy of P and P values, for those coun- 
tries which e.g. do not use barley, and thus their ordering 
becomes somewhat arbitrary). 

Let us discuss in more detail few concrete examples 
shown in Tables 1-5 in Appendix. For all commodities 
first 5 positions are very close in both ways of ranking. As 
a significant change we note Canada which moves from 
K* — 1\ down to K* — 16 and Mexico with respective 
change from K* — 13 to K* > 20: the export of these two 
countries is too strongly oriented on USA that becomes di- 
rectly visible through CheiRank analysis. In contrast Sin- 
gapore moves up from K* — 15 to K* — 11 that shows 
the stability and broadness of its export trade, a similar 
situation appears for India moving up from K* = 19 to 
K* = 12. 

Even more strong changes of ranking appear for spe- 
cific commodities. For example for crude petroleum Rus- 
sia moves up from K* = 2 to K* — 1 showing that its 
trade network in this product is better and broader than 
the one of Saudi Arabia. Iran moves in opposite direction 
from K* = 5 down to K* = 14 showing that its trade 
network is restricted to a small number of nearby coun- 
tries. A significant improvement of ranking takes place for 
Kazakhstan moving up from K* = 12 to K* = 2. The 
direct analysis shows that this happens due to an unusual 
fact that Kazakhstan is practically the only country which 
sells crude petroleum to the CheiRank leader in this prod- 
uct Russia. This puts Kazakhstan on the second position. 
It is clear that such direction of trade is more of politi- 
cal or geographical origin and is not based on economic 
reasons. 

For natural gas there are also significant differences 
between two ways of ranking. Thus, USA moves strongly 
up from K* — 10 to K* — 1 due its broad trade network in 
this product. Canada moves down from K* = 2 to K* — 8 
due to its too strong trade orientation on USA. A small 
country Trinidad and Tobago moves up from K* = 15 to 
K* = 2 since it provides about 70% of import of top leader 
USA. 

Significant reordering appears also for barley trade. 
Thus, the leader Ukraine moves down from K* = 1 to 
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Fig. 12. (Color online) Comparison of ranking between 
PageRank K and ImportRank K (left column), and between 
CheiRank K* and ExportRank K* (right column) for year 
2008. The shown commodities are; all commodities (panels a, 
b); crude petroleum (panels c, d); natural gas (panels e, f); 
barley (panels g, h); cars (panels i, j). Only top 100 ranks are 
shown. 



K* — 6 due to too narrow trade network and USA moves 
up from K* = 8 to K* = 3 due to its broad trade network. 

For trade of cars we have France going up from K* = 7 
to K* = 3 due to its broad export network. Also Thailand 
goes strongly up from K* = 19 to K* = 10 due to its 
broad trade links. We note that on the side of import we 
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Fig. 13. (Color online) Country positions in the ImportRank- 
ExportRank plane {K,K*) for year 2008. The shown com- 
modities are: all commodities (panels a, b); crude petroleum 
(panel c, d); natural gas (panel e, f); barley (panel g, h); cars 
in i) and j). Left column shows a global scale (227 countries) 
while right column illustrates the first 40 x 40 region. Data can 
be compared with those in Fig. [l] 



have strong change for Nigeria which moves from K > 20 
up to iiT = 1. This is the most populated country in Africa 
with a strong income due to oil trade which provides a 
large fraction of USA import. With such oil income Nige- 
ria buys cars from all over the world and thus becomes at 
the top of PageRank. 



Finally we note that among top 20 countries ranked 
in 2DRank K2 by all commodities in 2008 (see Table 1) 
there 14 among G-20 major economies [27 . At the same 
time ExportRank gives 13, and ImportRank gives 14 coun- 
tries from 19 of G-^O-list. We attribute a difference in 5 
countries to political and geographical aspects of G-20- 
selection. 



5 Discussion 

In this work we constructed the Google matrix of the 
WTN using the enormous UN COMTRADE database. 
From this matrix we obtained PageRank and CheiRank of 
all world countries in various types of trade products for 
years 1962 - 2009. This new approach gives a democratic 
type of ranking being independent of the trade amount of 
a given country. In this way rich and poor countries are 
treated on equal democratic grounds. In a certain sense 
PageRank probability for a given country is proportional 
to its rescaled import flows while CheiRank is proportional 
to its rescaled export flows inside of the WTN. 

The global characteristics of the world trade are an- 
alyzed on the basis of this new type of ranking. Even if 
all countries are treated now on equal democratic grounds 
still we find at the top rank the group of industrially de- 
veloped countries approximately corresponding to G-20 
(74%). Our studies establish the existence of two solid 
state domains of rich and poor countries which remain 
stable during the years of consideration. Other countries 
correspond to a gas phase with ranking strongly fluctuat- 
ing in time. We propose a simple random matrix model 
which well describes the statistical properties of rank dis- 
tribution for the WTN. 

The comparison between usual ImportRank-Export- 
Rank (see e.g. [3]) and our PageRank-CheiRank approach 
shows that the later highlights the trade flows in a new 
useful manner which is complementary to the usual analy- 
sis. The important difference between these two approaches 
is due to the fact that ImportRank-ExportRank method 
takes into account only global amount of money exchange 
between a country and the rest of the world while Page- 
Rank-CheiRank approach takes into account all links and 
money flows between all countries. We hope that this new 
approach based on the Google matrix will find further use- 
ful applications to investigation of various flows in trade 
and economy. 

Acknowledgments: We thank Gilles Saint-Paul (TSE, 
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Chepelianskaia (UN Bangkok) for insights on UN databases. 
Our special thanks are addressed to Arlene Adriano and 
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6 Appendix 



Table 1. Top 20 ranking for all commodities - 2008. 
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Table 2. Top 20 ranking for crude petroleum - 2008. 
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Table 3. Top 20 ranking for natural gas - 2008. 
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Table 4. Top 20 ranking for barley - 


2008. 
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Russian Fed. Russian Fed 


Israel 


Spain 


15 


Cyprus 


ludia 


Austria 


Colombia 


Bulgaria 


16 


Denmark 


Hungary 


Poland 


Algeria 


Netherlands 


17 


Occ. Palestinian Terr. Romania 


Brazil 


Kuwait 


Lithuania 


18 


Switzerland 


Belgium 


Ireland 


Brazil 


Sweden 


19 


Bosnia Herzegovina Lithuania 


France 


Morocco 


Belgium 


20 


Jordan 


Sweden 


South Africa Turkey 


India 
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Table 5. Top 20 ranking for cars 


- 2008. 




Ran 


K 


K* 




K 


K* 


1 


Nigeria 


GsriiicLiiy 




USA 


Germany 


2 


Germany 




USA 


Germany 


Japan 


3 


Prance 


USA 




UK 


USA 


4 


USA 




UK 


Prance 


Canada 


5 


Russian Fed. 




Belgiurn 


Italy 


Rep. of Korea 


6 


UK 


UK 


Spain 


Russian Ped. 


UK 


7 


Belgium 


Belgium 


Italy 


Belgium 


France 


8 


Ukraine 






Canada 


Spain 


9 


Italy 


kJ w LI- vi-X XV/ ci* 


A imtralia 


Spain 


Belgium 


10 


Greece 


Thailand 


Canada 


China 


Mexico 


11 


Venezuela 


Mexico 


China 


Netherlands 


Italy 


12 


Spain 


Italy 


Netherlands 


Australia 


Slovakia 


13 


China 


Canada 


U. Arab Emir. 


U. Arab Emir. 


Czech Rep. 


14 


Netherlands 


U. Arab Emir. 


Austria 


Saudi Arabia 


Poland 


15 


Australia 


Czech Rep. 


South Africa 


Austria 


Sweden 


16 


Japan 


Slovakia 


Poland 


Mexico 


Turkey 


17 


Albania 


Hungary 


Thailand 


Poland 


Hungary 


18 


Romania 


Australia 


Russian Ped. 


Switzerland 


Austria 


19 


Sudan 


Austria 


Turkey 


Finland 


Thailand 


20 


Canada 


China 


Portugal 


Ukraine 


South Africa 



